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We find model-independent upper limits on rates of dark matter annihilation in galactic halos. 
The Born approximation generally fails, while exotic threshold enhancements akin to "Sommerfeld 
factors" also turn out to be baseless. The most efficient annihilation mechanism involves perturba- 
tively small decay widths that have largely been ignored. Widths that are very small compared to 
TeV mass scales suffice to cause large enhancements in the velocity averaged cross sections. Bound 
state formation in weakly coupled theories produces small effects due to wave function normaliza- 
tions. Unitarity shows the Sommerfeld factor cannot produce large enhancements of cross sections, 
and serves to identify where those approximations break down. 
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I. THRESHOLD ENHANCEMENTS 

There is great interest in recent data from the 
PAMELA [J, FERMI and PPB-BETS @ experi- 
ments. The observations suggest a significant signal in 
excess positron production in galactic halos, as long sug- 
gested by the HEAT 4] and ATIC [| experiments. Pos- 
sible explanations range from exotic mechanisms @ , un- 
certain features of pulsars 0], to dark matter decays § 
and dark matter annihilation 

In considering annihilation there are puzzles from com- 
paring predictions of relic densities with rates of particle 
production in the current era. This has led to invoking 
more or less exotic threshold enhancements under the 
catch-phrase of "Sommerfeld factors" [l(| [H| • 

One reason to appeal to a Sommerfeld factor is to boost 
cross sections of TeV-scale particles well above Born- level 
estimates. We find the starting point of Born-level cross 
sections is not a good approximation for much different 
reasons. Basic facts of finite width particle physics sub- 
stantially revise estimates of annihilation rates in galac- 
tic halos. We find that annihilation of TeV-scale dark 
matter with typical electroweak couplings can actually 
saturate unitarity limits over the observable range. We 
also obtain upper limits to halo annihilation rates that 
do not depend on fine details of the dark matter velocity 
distribution. 

Non-relativistic scattering amplitudes can be classified 
by their analytic properties in the complex momentum 
plane. Stable bound states are described by poles on 
the positive imaginary axis. It follows that stable bound 
states produce no remarkable enhancement of annihila- 
tion rates in the physical region of real momentum k. 
Metastable particles or resonances, described by poles 
of finite width, are in no way comparable with stable 
bound states, because everything observable (and poten- 
tially large) is a strong function of the width. 

Unless one is considering an absolutely stable inter- 
mediate state, all intermediate states in particle physics 



have a finite lifetime. Pursuing the consequences of finite 
lifetimes with galactic halo kinematics re-directs atten- 
tion from exotic mechanisms to ordinary physics. There 
are two salient cases. If the width of an intermediate 
annihilation state is limited by the initial state veloc- 
ity v, then the peak of the cross section goes like 1/v 2 . 
This case produces the largest reaction rates in halos, 
and most conservative bounds. If the width of the inter- 
mediate state is constant, the peak of the cross section 
goes like 1/v. In these and intermediate cases the peak 
cross section actually dominates the entire halo velocity 
distribution for a surprisingly broad range of dark matter 
parameters. As a result, our more conservative bounds 
merge smoothly with reasonable estimates predicting sur- 
prisingly large rates. 



II. BREIT-WIGNER FORMULAS 

Relic particles trapped in galactic halos will be non- 
relativistic, with velocities v ~ 10 -3 . There are several 
distinctly different non-relativistic "Breit-Wigner " for- 
mulas. Most Breit-Wigner cross sections a res can be cast 
into the form 
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Here Bi and Bf are the branching fractions to the initial 
and final state, and k is the momentum of an initial state 
particle in the center of mass frame. Different values of 
the parameter N — 0, 1 distinguish two classic limits: 

Phase Space Limited Case, N = 0: It is common for 
2 -> 2 non-relativistic physics to be quasi-elastic. In 
particular, the final state phase space may be severely 
limited by the initial state velocity v. Ignoring spin and 
matrix elements, the Lorentz-invariant phase space in- 
tegral LLPS for two particles of momentum Pf,Pf and 
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Here vj is the final state velocity of either particle in 
the CM frame. When initial and final state masses are 
comparable, and the 2-body states dominate, the total 
width r ~ KVf ~ kv, where k absorbs coupling constants 
and matrix elements. Incorporating the explicit velocity 
dependence with an s-channel propagator leads to Eq. 
[T] with N — 0. Note that the peak of the cross section 
scales like 1 / (m 2 v 2 ) , making this case potentially capable 
of saturating elastic unitarity bounds. 

Relativistic Phase Space Case, N = 1: Anihilla- 
tion may also proceed to final states which are ultra- 
relativistic. Then the square root in Eq. [2] approaches 
1, and the partial width Tf ~ constant in this limit. 
Any other kinematic situation where Q 2 /mj goes to a 
finite constant as v — » will produce the same outcome. 
This includes the "exoergic" resonances long known in 
low-energy nuclear physics, and associated with the "1/v 
law" of low energy cross sections. These cross sections 
do not increase as fast as unitarity would allow as v — > 0. 

The difference between 1/v and 1/v 2 velocity depen- 
dence is dramatic. Yet is only part of the story, be- 
cause resonances may produce large cross sections ei- 
ther way. For example, neutron absorption cross sections 
on Gadolinium- 157 exceeding one hundred million barns 
have been observed, [l^ . This comes in the seemingly 
mild l/v case not impinging on a unitarity limit. The 
experimental stunt simply exploits neutrons with grossly 
small velocities of order 3 meters per second. In much 
the same way, galactic halo velocities of order 10~ 3 are 
grossly small on the scale of particle physics. The com- 
bination of low speed halo kinematics and very ordinary 
widths produces surprisingly large enhancements. 



III. GENERAL LIMITS ON THE VELOCITY 
AVERAGED BREIT-WIGNER CROSS SECTION 

The halo annihilation rate via a single s— wave reso- 
nance is governed by the velocity-weighted cross section 
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Here &haio(v) — dN/dv 3 is the normalized dark matter 
relative velocity distribution, assumed from astrophysics 
to be a smooth function on the scale of 100-500 km/s. 
In an isothermal halo model the velocity distribution is 
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While the actual velocity distribution is uncertain the 
phase space factors of v 2 are general. The isothermal 
halo will illustrate the method, but none of our upper 
bounds depend on it. 

The rate (o~ res v) res is a function of E , E res , T and 
mx- If other scales are expressed in units of m x ~ 
TeV the conjunction of several rapidly varying func- 
tions makes analysis troublesome, as noted by Griest and 
Seckel [l3j. However in the present universe the halo 
energy mxv 2 /2 ~ 10 _6 mx is rather small on particle 
physics scales. It is natural to rescale variables in units 
of the halo characteristic energy, defining 
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Assuming the equilibrium distribution, some algebra 
gives 



, , 2 2 - Ar (2 7 r)^V Ar - 1 , 

{aV)res = o ^ITO, ) , 

m x 



(4) 



where 
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Note that Zjv (7o> ^o) is analytic for all 70 > and 5q 
regardless of the sign of <5o- It can be computed exactly in 
terms of Exponential Integral (Ei) functions. We found 
it more useful to observe that 7jv(70) <^o) has certain ab- 
solute upper limits for all possible values of 70 > and 
S. Consider the derivative d Iq (70, So)/d^o- 
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Since the integrand above is positive definite, the integral 
achieves its maximum at 70 — > 00. For 70 >> 1 the 
integration becomes trivial, yielding Jjv(7oj So) < 1. A 
stronger limit notes the integrand of Eq. [5] is cut off 
for z < 7o when 70 < 1, S < I, implying 7jv(7o, S ) < 
1 — e~ c ~ <0 , where C is a constant. Numerical work shows 
that for all parameters 



/o(7o, *o) <l-e-*-» 
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These are close to equality for positive So << 1. Figure 
[T] shows a plot of In (lo, <5o)f° r a wide range of 70, So and 
how the integral approaches the upper bound. 
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FIG. 1: The integral Jjv(7o, So) (dark shaded) and upper 
limits cited in the text (transparent mesh). By Eq. [3]the rate 
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The positivity property of Eq. [6] holds for all halo 
distributions. The upper limit BW — > 1 produces a uni- 
versal inequality: 
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The expected value (1/v) is relative to the distribution 
<&haio( v )i n °t dN/dv. If the equilibrium distribution is 
assumed, then 
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The result is a possible significant enhancement fac- 
tor (EF) ("boost factor") for annihilation rates. The 
enhancement factor is denned relative to a typical Born 
approximation usorn — ^ct 2 x /m 2 x v 2 ~ N : 
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Note that the upper limit does not depend on the po- 
sition of the resonance nor on any halo properties. 



1. N = Enhancement Factors 

For N — 0, Eq. [10] leads to substantial enhance- 
ments approaching the unitarity bound when the fun- 
damental width r is large enough. Obtaining a "large 
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FIG. 2: Upper limits (diagonal lines) of resonantly enhanced 
annihilation rate (av) re3 in the isothermal halo distribution. 
Solid curves (black online) are computed with fixed r/mx. 
Gray triangle in upper right is the unitarity bound. The 
thick dashed curve (blue online) is the maximum value for 
the cross section for N = 1. Thin dashed curves (red on- 
line) show (<jv) res computed for bound state processes using 
F = a 5 x mx /2 and E res = — mxo? x /4,. Middle curve (or- 
ange online) is the neutrino-based upper limit of Ref. [17] . 
Horizontal line (green online) is a conventional lower bound 
(av) ~3x 10~ 26 cm 3 /s- 



enough" width from a weakly coupled theory might ap- 
pear special. Yet remember that halo annihilations are 
driven by the width in units of the rather small scale 
E ~ 10~ 6 mjf. For TeV-scale dark matter a width T > 
MeV is large enough to dominate the halo width and 
make BW(T, E res ) ~ 1. Recall that the J/ip has a width 
of order 0.1 MeV and is exceedingly "narrow". For an 
elementary particle on any mass scale of GeV-TeV not to 
have widths exceeding 10~ 6 mx requires special conspir- 
acies or selection rules. 

Figure [5] shows that even a tiny value of T/mx ~ 10~ 8 
can produce rates much larger than the oft-cited value 
(crv) ~ 3 x 10~ 26 cm 3 /s. It is a new insight that merely 
including physics of widths tends to saturate unitarity 
bounds in halo annihilation. 
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FIG. 3: Ratio of Ii/Iq (shaded area) compared to the uniform 
value of 1 (mesh) 



Yet just as above, everything about any significant 
enhancement depends strongly on the width, and won't 
proceed without it. To estimate widths, first note that 
bound states are spatially large for small coupling con- 
stant ax- The size of a weakly coupled bound state is 
roughly estimated by the "Bohr radius" do, where 

a ~ l/m x a X - 

Similarly, the binding energy is E res ~ mxa\. Next 
recall that the Schroedinger wave function at the origin 
■0(0) determines the width via T ~ \%jj(Q)\ 2 a c where <r c is 
a continuum cross section. 

The wave function at the origin is set by the inverse of 
the size of the bound state: 



\m\ 2 



2. N = 1 Enhancement Factors 

Equation [8] highlights a factor of (l/v) absent with a 
rclativistic phase space (N — 1). To a first approxima- 
tion the ratio of the N = 1 case relative to the N = 
case is of 0(vq). This is made more precise using Figure 
[3l which shows a plot of the calculated ratio of integrals 
I\ / Jo that remains. This ratio is of order unity for most of 
the parameter space, except the regions where 70 << 1. 
Once again, only when widths are very tiny do resonance 
widths not tend to swamp the halo distribution. 

While representing stronger limits, the bottom panel 
of Figure [2] again shows significant enhancements over 
a broad range of parameters of current interest. The 
difference between I\ and Jo tends to disappear whenever 
T/mx is not exceptionally small. In the next Section 
we turn to the metastable bound state case, which does 
happen to exhibit exceptionally small widths on general 
grounds. 



A. Metastable Bound States, and Narrow Resonances 

The case of annihilation passing through intermedi- 
ate metastable bound states has generated great interest. 
This case is different and deserves a separate discussion. 
Suppose dark matter interacts with a light messenger 
particle of mass /1, with coupling-squared ax- If the in- 
teraction is attractive, which is readily arranged for par- 
ticular spins, then non-relativistic physics predicts there 
is always a bound state for sufficiently large coupling. 
The conditions are 
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where k is a constant of order one. The demonstration is 
an easy variational calculation using a ground state Hy- 
drogenic wave function. A helpful discussion is also given 
in Ref. [l4j. For parameters m x ~ TeV and \i ~ GeV 
bound state formation needs ax > 10 -3 , which is well 
within the electroweak-scale couplings of most models. 



The continuum annihilation cross section a c ~ a 1 x ^ A for 
A > depending on the model. For reference the annihi- 
lation rates of ortho (para) positronium via three (two) 
photons go like a® m (a 5 em ). Thus bound state widths fol- 
low a general pattern 



m x 



< 10' 



The right hand side is a fair upper limit for ax ~ 10 ~ 2 . 
Restricted phase space factors and branching ratios can 
only reduce this. Comparing Eq ~ 10~ 6 mx, we find 
that 70 << 1 is by far the generic case for annihilation 
from a bound state. As a consistency check, consider 
the definite case of spin- 1/2 dark matter interacting with 
vector particles. Nature has already done this calculation 
with the J ftp decay via gluons, which has Tj/^/mj/^ < 
10~ 4 . The J/tp is sufficiently heavy that the perturbative 
phase space factors are driven by dimensional analysis, 
as expected for TeV-scale physics. The raw J/ip ratio 
needs to be re-scaled by (ax / a s ) 4+A ~ 10~ 4 , which gives 
satisfactory agreement. 

When 70 << 1 it is a good approximation to replace 
BW(T, E res ) ~ n(T/2)S(E - E res ). A short calculation 
then gives 
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where E, 



res — mxv^ es /2. This formula has no singularity 
as v res — > because &haio(vres) ~ Vres nas compensat- 
ing factors from phase space (Eq. [3]). If a metastable 
bound state resonance lies above threshold in an expected 
electroweak range the effects are quite small. Taking 
E res — mxa x /2 ~ 10~ 4 mx, and the equilibrium halo 
model with scale v = 10~ 3 , the factor e ~ Er,:s ^ Eo ~ e~ 100 
is too small to consider further. When the resonance is 
below threshold it must have width T > \E res \ to intrude 
into the physical region. Since T is proportional to sev- 
eral powers of ax compared to E res this case can also be 
set aside. If E res — ¥ with T >> E res is contemplated, 
it implies the decay time scale is much less than a binding 



5 



(orbital) time scale, which is not consistent with bound 
states forming in the first place. 

An exponentially small suppression can be avoided by 
adjusting the binding into the range probed by the halo 
velocity. For example choose ax ~ 10~ 3 . This device 
rapidly loses consistency because the bound state crite- 
rion ax ^ m/ to x needs couplings not too small. If a 
bound state is tuned to the vicinity of the peak, then 
the halo factors will be order unity. Meanwhile there re- 
mains in Eq. [TT] an overall factor of T/(mxce x ) << 1. 
Figure [3 compares the upper limits from annihilation of 
continuum processes (70 > 1 generically) to processes 
proceeding via the bound state (70 < 10 6 a^-) using the 
isothermal halo and conservative values BiBf — > 1. Vi- 
able enhancement mechanisms should also respect the 
neutrino-based bounds of Mack, Beacom and Bell [I?} 
included in the Figure. In case of N = 0, a bound state 
could cause large cross section enhancements, but only 
for couplings ax > 0.1 which are beyond the stable per- 
turbative regime. In case of N = 1 the limits for bound 
states are even tighter. 

Figure [2] shows that a single bound state with per- 
turbative couplings has no chance of causing significant 
enhancements. Except for strong coupling, there is no 
dynamical mechanism to generate large enhancement fac- 
tors from non-relativistic bound state resonances in the 
current universe. The conclusion does not depend on 
the spin or quantum numbers of new physics, and is too 
strong to escape by adding up several resonances, unless 
they are so numerous their numbers alone overcome small 
couplings, as for KK modes. 



1. Breit Wigner Effects on Relic Abundance 

Relic abundance is a different topic than halo annihi- 
lation. Ibe, Murayama, and Yanagida 15[,and Guo and 
Wu [l6| have calculated thermal evolution for the case 
of a narrow state close to threshold. Their model cross 
section is essentially equivalent to our N — 1. The reso- 
nance position is close enough to threshold for its effects 
to overlap into the physical region during relic evolution. 
They find that even a tiny ratio of width to resonance 
invariant mass, denoted by 7 = T/M res < 10 -3 , pro- 



duces significant effects on relic densities compared to 
traditional constant cross section approximations. The 
sense of this effect causes a relative decrease in annihila- 
tion rates in the early universe, which tends to leave too 
much relic. To keep the relic density Qx fixed, avoid- 
ing over-closure, Refs. [H[ introduce "boost factors" 
to correct the normalization parameters of the cross sec- 
tion o~o ~ BiBf/m x . When those boost factors are ap- 
plied directly to halo annihilation, they develop much 
larger cross sections than the well-known cosmological 
value erg ~ 10~ 9 GeV~ 2 , which might be relevant to the 
PAMELA-ATIC observations. 

However, in general it is not possible to go from the 
relic calculation to the halo calculation directly in this 



manner. The halo annihilation rate < o~v > has a new 
and separate sensitivity that is a priori disconnected from 
relic calculations. The halo estimates are driven by the 
new parameter 70 = j M res /2Eo. As long as 70 >> 1, 
the upper bounds on the halo annihilation cross section 
will be saturated. This key feature is conceptually absent 
if the halo annihilation cross section is simply re-scaled 
by factors invoked for relic evolution. Thus the reported 
"boost factors" of the relic calculations do not take into 
account the Breit- Wigner effects on halo annihilation we 
have found. This explains why the suggestion [TH, [iH ] 
that very small T/M res is necessary or tends to enhance 
halo rates is not general, and appears different from our 
conclusion. It is clearly possible to find models and pa- 
rameter regions where both, or neither of the correct relic 
density and halo enhancement phenomena can be accom- 
modated. 

We note that relic densities are also subject to many 
uncertainties of galaxy formation and the other boost 
factors representing "dumpiness" . For purposes of con- 
fronting experimental data, it seems best to separate the 
problems of halo annihilation and relic evolution entirely, 
despite mathematical similarities in how they are calcu- 
lated. 



IV. SOMMERFELD FACTORS 

We have shown that Breit- Wigner width effects of typ- 
ical particle physics type can be surprisingly large, while 
bound state effects have little chance to compete. Som- 
merfeld factors have also been claimed as a mechanism to 
produce large enhancements not involving particle widths 

Given an s— wave cross section er°, which has been 
computed in the plane wave basis, the Sommcrfcld-based 
recipe to include Coulomb wave effects is to make a re- 
placement 



a°S(v, a); 



S(v, a) = — 



2ir 



v \ _ q—2-ko./v ' 

Here a is the fine structure constant. Since cross sections 
contain many other terms of different orders in a? /v k , to- 
gether with logarithmic type dependence, the recipe is an 
approximation by re-summation of selected contributions 

M. 



2. Motivation for Re- summation 

Non-relativistic QED has complicated logarithmic and 
power-behaved infrared singularities. Singular terms 
must be summed or controlled in some way to avoid up- 
setting perturbation theory. There are reasons to believe 
that the leading singularities appear [28| order by order 
as a series in a J /v J . Evidently such a series is summed 
in S(v, a) = S(a/v). 
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Sub-leading terms are dropped in any re-summation - 
for example, a term of order a J /v k with j > k is sub- 
leading as v — > 0. The fact that infinitely many sub- 
leading terms exist comes from the fact that it is always 
possible to add a photon exchange loop to any diagram, 
and all loops have some integration region not singular 
as v — > 0. 

The purpose of re-summation is to extend the reach of 
perturbation theory into the difficult, non-perturbative 
regime. Some examples illustrate typical limitations. Ex- 
pand S(v, a) ~ l+ir(a/v)+n 2 a 2 /3v 2 +.... For a = 1(T 2 , 
and v = 5 x 1(T 2 the third term is tt 2 /75 ~ 0.13. It 
happens to be larger than a typical non-singular term 
of order a; retaining it is well-motivated for these kine- 
matics. The series is also stable in this regime: a high 
order term (ira/v) 10 ~ 0.0096 is small. Yet at the 
smaller velocity v = 5/1000, the sub-leading correction 
a 11 (ir/v) 10 — > 958956.0 is not included, is hardly small, 
and the leading-order re-summation fails. It is the nature 
of such re-summation that self-consistency breaks down 
in the region a/v > 1, exactly where S(v, a) > 1. 

If one believes a re-summation recipe might generate 
corrections of order 30-50%, say, there's seldom any rea- 
son to invoke it for factors of "10" or more. Yet recent 
treatments of dark matter annihilation have imagined the 
Sommerfeld factor to be very general. It has been held 
responsible for extremely large enhancement factors of 
S » 10, while also coming from any generic interac- 
tion involving light Yukawa particles [20( . The perception 
comes from a practice of citing continuum Coulomb nor- 
malization factors |-0c 1 2 (0) . Since Coulomb normaliza- 
tions are known exactly, the procedure has been thought 
to be "exact." 

We have traced early literature to find several logi- 
cal and historical contradictions. Guth and Mullin high- 
lighted the approximations in 1951 [21] . In lowest order 
approximation, but while using Coulomb wave functions 
ipci for a basis, one encounters matrix elements M of the 
form 

M= J d 3 xr C 2V^a- 

Insert complete sets of momentum eigenstates jfc): 

M = J d 3 kd 3 k' ^* C2 {k)V kk ^ C i{k'). 

The Coulomb wave functions are sharply peaked in the 
vicinity of certain momenta k\, fc2, identified by taking 
the limit a — 0. Assume the plane-wave matrix elements 
Vk lt fc 2 = (ki I ^1^2) are relatively smooth functions of mo- 
mentum transfer. Make a rough approximation moving 
^fcifc^ outside the integral: 

M ->• V kuk2 J d 3 kr C 2(k) J d 3 k'vj ci {k), 

V kuk2 iP* C 2(x = 0)tPc i(x = 0). (12) 

In the last line ipc( x — 0) appears as the coordinate- 
space wave function at the origin, "improving" the plane 



wave calculation. Inserting the analytically known nor- 
malization of \ipc\ 2 then produces the factor S(v, a) for 
the cross section. 

The operations of separating the collision and wave 
function integration into products is one of leading power 
factorization. It is used in QCD calculations separating 
"hard" and "soft" regions of perturbation theory, but 
in that case while attempting to be systematic. Care- 
ful work with positronium annihilation [22[ does not use 
the factorized approximation. Instead, reference to re- 
summation is made after the full calculations are car- 
ried out. An early work [22[ on one-loop corrections to 
positronium decay states that "Coulomb effects are in- 
cluded by this (factored) method to all orders in e 2 , 
though only, of course, approximately." (Italics are ours.) 

What did Sommerfeld actually do? We consulted 
his 1931 article, in German, to see it introduced ex- 
act Coulomb wave functions to calculate bremsstrahlung, 
while it never suggested factorization. It is a tour de force 
of early quantum theory; consulting it for a renormaliza- 
tion factor actually perpetuates a normalization mistake. 
Cross sections are defined by ratios relative to a flux com- 
puted with a given normalization. The overall normaliza- 
tion of physical states cancels out in total cross sections: 
and so Eq. [T2]is not only approximate, it is incomplete. 

Elwert's 1939 dissertation [23| recognized this, as 
as did Guth in 1941 24]. These papers abandoned 
Sommerfeld's calculation and used the ratio of two in- 
and out- Coulomb factors as an approximate factorized 
ansatz. "Elwert factors" are used in atomic and molec- 
ular physics to cancel spurious pre-factors going like v 
from other approximations, but only when their effects 
are not too large. Experimental confrontation of the Elw- 
ert factor finds errors of relative order unity in the region 
the factors are of order unity [H, [25| . Elwert and collab- 
orators find this kind of breakdown reasonable (23|. In 
no event are very large corrections ever credited. 



3. Multiplicative Factors Must Fail 

In retrospect, we find the concept behind generating 
singularities via multiplicative factors questionable on 
general grounds. 

A general scattering amplitude has the partial wave 
expansion 

m fc ) = \ Y,^ 2l + l )^ k ) p i { - cos9 )- 
1 

For each partial wave cross section o~i of angular momen- 
tum /, elastic unitary gives the upper limit 

47r(2Z + l)|/ ; (fc)| 2 4tt(2Z + 1) 
° l k 2 k 2 

This summarizes the unitarity bound of Ref. [26| . Since 
each partial wave has a finite cross section, no partial 
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wave can possibly have a singularity. "Improving" the s- 
wave cross section - or any particular partial wave cross 
section - by singular terms of order ax / v then contradicts 
unitarity for ax/v ^ 1< This is just the same region 
where the claimed Sommerfeld factor S(v, a) >> 1. 

Can one escape the contradiction by appealing to small 
ax? It seems not: No small value for ax is used in the 
logic of an exact normalization citing Sommerfeld. Small 
coupling is also no protection from internal inconsistency. 
Unitarity and analyticity in perturbation theory are ex- 
act facts maintained in a systematic way, order by order, 
regardless of the size of the coupling constant, small or 
large. When violated, it shows the calculation was bad, 
just as indicated by sub- leading terms [20j| . 

This problem of consistency is different from the 
one previously recognized. Dark matter interactions 
have finite range, while the infrared singularities of re- 
summation come from infinite range. To account for 
this the authors of Ref. [ll| argued that for a finite 
range potential, a Sommerfeld enhancement would satu- 
rate when the deBroigle wavelength of colliding particles 
would be larger than the range of the force. A related 
statement actually follows from a WKB approximation: 
when the de Broglie wavelength is tiny compared to the 
range, and one works with wave functions at short dis- 
tance, the range effects of a Yukawa potential e _A ""/r 
drop out. Note that the range criteria do not depend 
on a coupling constant, and also don't specify any par- 
ticular angular momentum channels. Yet the singulari- 
ties of scattering amplitudes, and particularly Coulomb 
singularities, do depend on the couplings and angular 
momentum channels. Whatever the scale where analo- 
gies between massless and massive models break down, 
the facts of partial wave unitarity are more general, and 
take precedence. They preclude large enhancement in 
any particular channel of fixed angular momentum. 

This leads to another useful bound. Replace every par- 
tial wave by the one with largest magnitude \f ma x\- Sum 
them up: The result is 

4tt^ 2 ^\f max \ 2 (l max + l) 2 

a l^y 2l + l )\j™x\ = p ( 13 ) 

This is a strict upper limit. The Sommerfeld factor is not 
a resonance, and in the absence of resonances \f m ax\ 2 
for every partial wave is small in weakly coupled theo- 
ries, making a small. The notion of canceling a pertur- 
bative factor of a\ with an enhancement of 10 4 (say) 
needs l ma x >> 10 2 partial waves. This is difficult to 
conceive - or at least a high burden of proof - for finite- 
range, perturbatively coupled Yukawa models. It is even 
more problematic that the angular momentum involved 



in annihilation is strictly limited by the quantum num- 
bers of intermediate states. When the intermediate state 
consists of a single particle, the angular momentum is 
bounded by its spin, closing the door on large spin sums. 

Since each partial wave is finite, how do some 
Coulomb-dominated processes, such as Rutherford scat- 
tering, actually become singular as v — > 0? As Wigner 
[27j and many others have noted, the Coulomb singular- 
ity is very special. On semi-classical grounds (actually 
the facts of Legendre series), in Eq. [13] the upper limit 
lmax ~ r max k, where r max is the range of the potential. 
This gives 

a < Airr 2 If I 2 < Airr 2 

The Coulomb singularity occurs because (1) the effec- 
tive range r max — > 00, and (2) an infinite number of 
partial waves actually can contribute. Closely related is 
Wigner 's classic theorem [27jj that power-law potentials 
V(r) > 1/r 2 are needed to develop any kind of singular- 
ity. 



V. CONCLUDING REMARKS 

We have explored significant effects in halo annihila- 
tion rates due to natural widths of intermediate states. 
The problem is intricate due to subtle interplay of energy 
scales. The Born approximation almost always fails, giv- 
ing gross underestimates of reaction rates. Tiny values of 
galactic halo velocities reverse an assumption that prop- 
agator widths might be "small corrections." That per- 
ception comes from comparing widths to particle masses, 
and does not capture the important features of halo an- 
nihilation. 

Given that TeV-scale particles with typical electroweak 
type couplings may easily have T/E » 1, Breit- Wigner 
factors of ordinary radiative corrections must generally 
be taken into account. Consistency of rates in particular 
channels, such as the apparent dominance of leptons, still 
needs to be considered model by model. The fact that 
merely including basic physics of widths may be quite 
significant. It revises the basic picture of annihilation 
scenarios confronting the PAMELA-FERMI-PPB-BETS 
data in a positive way that increases the possibilities to 
find new physics. 
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